Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?

Seascape Genomics and Phylogeography of the Sailfish (Istiophorus platypterus).

Permeable phylogeographic barriers characterize the vast open ocean, boosting gene flow and counteracting population differentiation and speciation of widely distributed and migratory species. However, many widely distributed species consists of distinct populations throughout their distribution, evidencing that our understanding of how the marine environment triggers population and species divergence are insufficient. The sailfish is a circumtropical and highly migratory billfish that inhabits warm and productive areas. Despite its ecological and socioeconomic importance as a predator and fishery resource, the species is threatened by overfishing, requiring innovative approaches to improve their management and conservation status. Thus, we presented a novel high-quality reference genome for the species and applied a seascape genomics approach to understand how marine environmental features may promote local adaptation and how it affects gene flow between populations. We delimit two populations between the Atlantic and Indo-Western Pacific oceans and detect outlier loci correlated with sea surface temperature, salinity, oxygen, and chlorophyll concentrations. However, the most significant explanatory factor that explains the differences between populations was isolation by distance. Despite recent population drops, the sailfish populations are not inbred. For billfishes in general, genome-wide heterozygosity was found to be relatively low compared to other marine fishes, evidencing the need to counteract overfishing effects. In addition, in a climate change scenario, management agencies must implement state-of-the-art sequencing methods, consider our findings in their management plans, and monitor genome-wide heterozygosity over time to improve sustainable fisheries and the long-term viability of its populations.

Genomics reveals the role of admixture in the evolution of structure among sperm whale populations within the Mediterranean Sea.

In oceanic ecosystems, the nature of barriers to gene flow and the processes by which populations may become isolated are different from the terrestrial environment, and less well understood. In this study we investigate a highly mobile species (the sperm whale, Physeter macrocephalus) that is genetically differentiated between an open North Atlantic population and the populations in the Mediterranean Sea. We apply high-resolution single nucleotide polymorphism (SNP) analysis to study the nature of barriers to gene flow in this system, assessing the putative boundary into the Mediterranean (Strait of Gibraltar and Alboran Sea region), and including novel analyses on structuring among sperm whale populations within the Mediterranean basin. Our data support a recent founding of the Mediterranean population, around the time of the last glacial maximum, and show concerted historical demographic profiles in both the Atlantic and the Mediterranean. In each region there is evidence for a population decline around the time of the founder event. The largest decline was seen within the Mediterranean Sea where effective population size is substantially lower (especially in the eastern basin). While differentiation is strongest at the Atlantic/Mediterranean boundary, there is also weaker but significant differentiation between the eastern and western basins of the Mediterranean Sea. We propose, however, that the mechanisms are different. While post-founding gene flow was reduced between the Mediterranean and Atlantic populations, within the Mediterranean an important factor differentiating the basins is probably a greater degree of admixture between the western basin and the North Atlantic and some level of isolation between the western and eastern Mediterranean basins. Subdivision within the Mediterranean Sea exacerbates conservation concerns and will require consideration of what distinct impacts may affect populations in the two basins.

Range-wide whole-genome resequencing of the brown bear reveals drivers of intraspecies divergence.

Population-genomic studies can shed new light on the effect of past demographic processes on contemporary population structure. We reassessed phylogeographical patterns of a classic model species of postglacial recolonisation, the brown bear (Ursus arctos), using a range-wide resequencing dataset of 128 nuclear genomes. In sharp contrast to the erratic geographical distribution of mtDNA and Y-chromosomal haplotypes, autosomal and X-chromosomal multi-locus datasets indicate that brown bear population structure is largely explained by recent population connectivity. Multispecies coalescent based analyses reveal cases where mtDNA haplotype sharing between distant populations, such as between Iberian and southern Scandinavian bears, likely results from incomplete lineage sorting, not from ancestral population structure (i.e., postglacial recolonisation). However, we also argue, using forward-in-time simulations, that gene flow and recombination can rapidly erase genomic evidence of former population structure (such as an ancestral population in Beringia), while this signal is retained by Y-chromosomal and mtDNA, albeit likely distorted. We further suggest that if gene flow is male-mediated, the information loss proceeds faster in autosomes than in X chromosomes. Our findings emphasise that contemporary autosomal genetic structure may reflect recent population dynamics rather than postglacial recolonisation routes, which could contribute to mtDNA and Y-chromosomal discordances.

A chromosome-scale reference genome assembly of the great sand eel, Hyperoplus lanceolatus.

Despite increasing sequencing efforts, numerous fish families still lack a reference genome, which complicates genetic research. One such understudied family is the sand lances (Ammodytidae, literally: "sand burrower"), a globally distributed clade of over 30 fish species that tend to avoid tidal currents by burrowing into the sand. Here, we present the first annotated chromosome-level genome assembly of the great sand eel (Hyperoplus lanceolatus). The genome assembly was generated using Oxford Nanopore Technologies long sequencing reads and Illumina short reads for polishing. The final assembly has a total length of 808.5 Mbp, of which 97.1% were anchored into 24 chromosome-scale scaffolds using proximity-ligation scaffolding. It is highly contiguous with a scaffold and contig N50 of 33.7 and 31.3 Mbp, respectively, and has a BUSCO completeness score of 96.9%. The presented genome assembly is a valuable resource for future studies of sand lances, as this family is of great ecological and commercial importance and may also contribute to studies aiming to resolve the suprafamiliar taxonomy of bony fishes.

Detecting genetic signals of selection in heavily bottlenecked reindeer populations by comparing parallel founder events.

Founder populations are of special interest to both evolutionary and conservation biologists, but the detection of genetic signals of selection in these populations is challenging due to their demographic history. Geographically separated founder populations likely to have been subjected to similar selection pressures provide an ideal but rare opportunity to overcome these challenges. Here we take advantage of such a situation generated when small, isolated founder populations of reindeer were established on the island of South Georgia, and using this system we look for empirical evidence of selection overcoming strong genetic drift. We generated a 70 k ddRADseq single nucleotide polymorphism database for the two parallel reindeer founder populations and screened for signatures of soft sweeps. We find evidence for a genomic region under selection shared among the two populations, and support our findings with Wright-Fisher model simulations to assess the power and specificity of interpopulation selection scans-namely Bayescan, OutFLANK, PCadapt and a newly developed scan called Genome Wide Differentiation Scan (GWDS)-in the context of pairwise source-founder comparisons. Our simulations indicate that loci under selection in small founder populations are most probably detected by GWDS, and strengthen the hypothesis that the outlier region represents a true locus under selection. We explore possible, relevant functional roles for genes in linkage with the detected outlier loci.

SambaR: An R package for fast, easy and reproducible population-genetic analyses of biallelic SNP data sets.

SNP data sets can be used to infer a wealth of information about natural populations, including information about their structure, genetic diversity, and the presence of loci under selection. However, SNP data analysis can be a time-consuming and challenging process, not in the least because at present many different software packages are needed to execute and depict the wide variety of mainstream population-genetic analyses. Here, we present SambaR, an integrative and user-friendly R package which automates and simplifies quality control and population-genetic analyses of biallelic SNP data sets. SambaR allows users to perform mainstream population-genetic analyses and to generate a wide variety of ready to publish graphs with a minimum number of commands (less than 10). These wrapper commands call functions of existing packages (including adegenet, ape, LEA, poppr, pcadapt and StAMPP) as well as new tools uniquely implemented in SambaR. We tested SambaR on online available SNP data sets and found that SambaR can process data sets of over 100,000 SNPs and hundreds of individuals within hours, given sufficient computing power. Newly developed tools implemented in SambaR facilitate optimization of filter settings, objective interpretation of ordination analyses, enhance comparability of diversity estimates from reduced representation library SNP data sets, and generate reduced SNP panels and structure-like plots with Bayesian population assignment probabilities. SambaR facilitates rapid population genetic analyses on biallelic SNP data sets by removing three major time sinks: file handling, software learning, and data plotting. In addition, SambaR provides a convenient platform for SNP data storage and management, as well as several new utilities, including guidance in setting appropriate data filters. The SambaR source script, manual and example data set are distributed through GitHub: https://github.com/mennodejong1986/SambaR.

Demography and adaptation promoting evolutionary transitions in a mammalian genus that diversified during the Pleistocene.

Species that evolved in temperate regions during the Pleistocene experienced periods of extreme climatic transitions. Consequent population fragmentation and dynamics had the potential to generate small, isolated populations where the influence of genetic drift would be expected to be strong. We use comparative genomics to assess the evolutionary influence of historical demographics and natural selection through a series of transitions associated with the formation of the genus Capreolus, speciation within this genus during the Quaternary and during divergence among European roe deer (C. capreolus) populations. Our analyses were facilitated by the generation of a new high-coverage reference genome for the Siberian roe deer (C. pygargus). We find progressive reductions in effective population size (Ne ), despite very large census sizes in modern C. capreolus populations and show that low Ne has impacted the C. capreolus genome, reducing diversity and increasing linkage disequilibrium. Even so, we find evidence for natural selection shared among C. capreolus populations, including a historically documented founder population that has been through a severe bottleneck. During each phylogenetic transition there is evidence for selection (from dN/dS and nucleotide diversity tests), including at loci associated with diapause (delayed embryonic development), a phenotype restricted to this genus among the even-toed ungulates. Together these data allow us to assess expectations for the origin and diversification of a mammalian genus during a period of extreme environmental change.